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Abstract 

In chaotic dynamical systems, an infinitesimal perturbation is exponentially 
amplified at a time-rate given by the inverse of the maximum Lyapunov ex- 
ponent A. In fully developed turbulence, A grows as a power of the Reynolds 
number. This result could seem in contrast with phenomenological arguments 
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suggesting that, as a consequence of 'physical' perturbations, the predictabil- 
ity time is roughly given by the characteristic life-time of the large scale struc- 
tures, and hence independent of the Reynolds number. We show that such 
a situation is present in generic systems with many degrees of freedom, since 
the growth of a non-infinitesimal perturbation is determined by cumulative 
effects of many different characteristic times and is unrelated to the maximum 
Lyapunov exponent. Our results are illustrated in a chain of coupled maps 

and in a shell model for the energy cascade in turbulence. 
05.45.+b,47.27.Gs 
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I. INTRODUCTION 

After the seminal work of Lorenz 0, it is well understood that the predictability of the 
state of a system ruled by a deterministic evolution law has severe limitations in presence of 
deterministic chaos. In systems with sensitive dependence on the initial condition one has 
an exponential divergence of the distance 5x between two initially close trajectories, i.e. 

\6x(t)\ ~ \Sx(0)\e M (1) 

where A is the maximum Lyapunov exponent [H. Consequently, if |5x(0)| = Sq and one 
accepts a maximum tolerance S max on the knowledge of the state of the system, (|l|) implies 
that the systems is predictable up to a time 

r (2) 

Equation (|2]) gives only a first rough answer to the problem since it does not take into 
account some important features of chaotic systems. Endeed to study the predictability of 
a generic dynamical system one has to consider the following non trivial aspects: 

• a) The Lyapunov exponent A is a global quantity: it measures the average exponential 
rate of divergence of nearby trajectories. In general there are finite-time fluctuations of 
this rate and it is possible to define an 'instantaneous' rate 7, called effective Lyapunov 
exponent ||, which depends on the particular point of the trajectory x(i) where the 
perturbation is performed. In the same way, the predictability time T fluctuates, 
following the 7-variations. 

• b) In dynamical systems with many degrees of freedom, the interactions among differ- 
ent parts of the system play an important role on the growth of a perturbation. The 
statistics of the effective Lyapunov exponent is not sufficient to characterize the growth 
of infinitesimal perturbations and one has to analyze the behavior of the tangent vec- 
tor z(t), i.e. of the direction along which an infinitesimal perturbation grows, see e.g 
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[|]]. Moreover, if one is interested in the behaviour of a perturbation concentrated on 
certain degrees of freedom, e.g. small length scales in weather forecasting, and in a 
prediction on the evolution of other degrees of freedom, e.g. large length scales, a 
relevant quantity is the time Tr necessary to the tangent vector to relax on the time 
dependent eigenvector e(t) of the stability matrix, corresponding to the maximum 
Lyapunov exponent. If the perturbations are small enough, i.e. £x is proportional to 
the tangent vector, one has that 

r~r, + ib,*= (3) 

A 6q 

where in general Tr may depend on <5x. So the mechanism of transfer of the error 5x 
through the degrees of freedom of the system could be more important than the rate 
of divergence of nearby trajectories ||. 

• c) In systems with many characteristic times, such as the eddy turn-over times in fully 
developed turbulence, if the perturbation is not infinitesimal, or if the threshold of 
accepted error is not small, T is determined by the detailed mechanism of transfer due 
to the non-linear effects in the evolution equation for 5x. In this case, the predictability 
time might have no relation with the Lyapunov exponent and T depends in a non- 
trivial way on the details of the system. 

The aspects a) and b) have been studied in previous works [0-0. In this paper we 
mainly address point c). We investigate the predictability problem for non-infinitesimal 
perturbations in the framework of two models: a system of coupled maps, and a shell model 
for the energy cascade in three dimensional turbulence. 

The paper is organized as follows. In section [TI| we discuss some phenomenological results 
in fully developed turbulence. In section [Til] the predictability problem is discussed in a 
system of coupled maps, using analytic methods supported by some numerical investigations. 
Section [TV] analyzes the more realistic case of a shell model of turbulence by both numerical 
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simulations and closure approximations. In section [V| we discuss the physical relevance of 
our results and open problems. 



II. SOME RESULTS FOR THE PREDICTABILITY IN FULLY DEVELOPED 

TURBULENCE 

In three dimensional fully developed turbulence, the inverse maximum Lyapunov expo- 
nent is roughly proportional to the smallest characteristic time of the system, the turn-over 
time r of eddies of the size of the Kolmogorov length rj. The argument, due to Ruelle ||, is 
the following. The Kolmogorov theory predicts that the longitudinal velocity difference at 
distance I = |r| scales as 

/ / \ 1/3 

v{£) = | v(x + r) - v(x) | ~ e 1 / 3 ^ 3 ~V ij-j (4) 

where V and L are the typical velocity and length of the energy containing eddies, and e is 
the mean rate of energy dissipation. 

The non-linear transfer of energy is stopped at the Kolmogorov scale rj where viscosity 
v is able to compete with the convective term, i.e. v ~ r)v(r]), thus from M) we have 



V 



LRe-V 4 (5) 



where Re = VL/v is the Reynolds number. The corresponding turn-over time is 

where To = L/V is the life time of the large scale disturbances. 

These dimensional relations imply that the maximum Lyapunov exponent scales with 
Re as 

A ~ — ~ — Re 1/2 (7) 

Taking into account the intermittency one expects that the presence of quiescent quasi- 
laminar periods changes the chaotic features of the fluid flow. The intermittency of energy 



dissipation can be described by introducing a spectrum of singularities h of the velocity 
gradients ||, i.e. assuming a local scaling invariance so that in any point x of the fluid 
v(£) ~ i h w . In the framework of the multifractal approach, one thus finds R 



A ~ — Re a with a = max 
T h 



D(h) -2-h 



(8) 



l + h 

where D(h) is the fractal dimension of the set of fluid points characterized by a given 
singularity h. The value of a depends on D(h). By using the function D(h) obtained by 



fitting the exponents (, q of the velocity structure functions with the random beta model |fL0 
one has a = 0.459..., slightly smaller than the Ruelle prediction a = 0.5, see ([7]) 

Relation (|2|) and (|8|) tell us that considering a very small perturbation at t — and a very 
small tolerance 5 max , the predictability time vanishes as the Reynolds number increases, as 
verified in recent numerical simulations of the shell model (see sect. [TVp ||. 

On a very different ground, without considering the exponential growth of infinitesimal 
disturbances, Lorenz [0, see also fl2| and |^3| , proposed a phenomenological approach to 



the predictability problem in turbulence. Consider wave-numbers around k with correspond- 
ing typical spatial scale I ~ k" 1 . The time r(k) for a perturbation at wavenumber 2k to 
induce a complete uncertainty on the velocity field on the wave number k is assumed to be 
proportional to the typical eddy turn-over time at scale £, from 

£ ( £ \ 2/3 

An incertitude 0{v{rj)) propagates through an inverse cascade from the Kolmogorov scale 
r] ~ k^ 1 up to the scale of the energy containing eddies L ~ . The predictabilty time on 
the large scales is therefore 

N 

t = Y, T(2 n v ) (io) 

ra=0 

where iV = log 2 (/cd/fco) ~ log Re. The geometric series ([TOf) is dominated by the term n = N 
so that T is practically independent of Reynolds number: 

T ~ T ~ L/V. (11) 



It is worth remarking that closure approximations |12|]T3|j allows one to write down simple 
equations for the evolution of disturbances in turbulent flows and fully confirm these results. 
In appendix B, we derive this type of equations in the simplified framework of the shell 
model, but there are no conceptual differences with the analogous derivation in the case of 
the Navier-Stokes equations. 

The Lorenz picture of an inverse cascade of a perturbation also permits to estimate 
the growth of a perturbation at intermediate times. Indeed, after a time t, a perturbation 
localized on the Kolmogorov length scale, is expected to affect an eddy of size It such that 
its turn-over time is r(£t) ~ t. It follows lfl4| , see @ 



t \ 3 / 2 

L ( 4- ) (12) 



.To, 

Let us indicate by v the reference velocity field and by v' the perturbed field. The size of 
the difference 5v(t) between two velocity fields at time t can be estimated as the typical 
velocity of a disturbance at scale £ t , so that from (13) and (O) 



\5v(t)\~v(e t )~V^t/T . (13) 

Perturbations grow as square root of time. Let us stress again that such perturbations 
cannot be considered infinitesimal so that the inverse cascade picture is not in contrast 
with the exponential amplification of errors on initial conditions that is present in chaotic 
dynamical systems. 

There are two main predictions of the Lorenz approach: predictability time on large scale 
independent of Reynolds number and growth of perturbations as a power of time instead of 
an exponential. 

In the next sections we use the theory of dynamical systems with many degrees of freedom 
to argue that while the former is correct, the latter is a too rough description of the real 
behavior of 5v(t). 
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III. PREDICTABILITY IN A CHAIN OF COUPLED MAPS 



In this section we discuss the predictability problem in a system of coupled maps each 
with a different time scale. This is an idealized situation where each degree of freedom has 
a chaotic dynamics with its own characteristic time and it is coupled to the other degrees 
of freedom by a weak local interaction. Our toy model can be considered the prototype 
for physical situations where one can separate the evolution on different scales. Despite its 
simplicity, it displays non-trivial properties which enlighten the behaviour of more complex 
systems. 

The system is given by a chain of chaotic coupled maps 

Xl (t + 1) = (1 - e)g l {x l (t)) + ef(x 2 (t)) 
x 2 (t + 1) = (1 - e)g 2 (x 2 (t)) + ef(x 3 (t)) 

(14) 

x N _i{t + 1) = (1 - e)g N _ 1 (x N „ 1 (t)) + ef(x N (t)) 
x N (t + 1) = g N (x N (t)) 

where the N variables Xk are defined on the interval J = [0, 1], e is the coupling constant, 
gk and / are functions of I into itself. In order to mimic the energy cascade of turbulent 
flow, we assume that the function g^ represents the evolution law of an eddy of size 4 in the 
absence of interactions. In this context, x k can be regarded as the velocity difference of the 
eddies on the octave of length scale = L 2~ k . The evolution is chaotic with a Lyapunov 
exponent proportional to the inverse of the eddy turn-over time r^. The simplest possible 
choice for gk{x) is the piecewise linear map 

9k(x) = e Xk x mod 1 (15) 

From dimensional arguments, and using (§) we can estimate the eddy turn-over time r k 
so that the Lyapunov exponent is 

\ k = ±2 2 ^ (16) 
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The coupling parameter e is small and constant for each scale, although more realistic 
models could consider scale-dependent couplings. The form of f(x) - the term representing 
the local interaction between eddies - does not affect the qualitative results. We have 
considered two possibilities: f(x) = x and f(x) = x 2 . 

In order to study the predictability, the system (|TJp has been integrated with two different 
initial conditions x(0) and x'(0) differing on the smallest scale, A&(0) = x' k (0) — Xk{0) = 
5k } N^o- The initial incertitude is thus confined to the smallest and fastest scale xn while we 
want to forecast the behavior of the system at the largest scale x\. 

Because time-scales are well separated, we expect that at the beginning the uncertainty 
in the system is driven by the fastest scale x^ and the error at any scale grows with the 
smallest characteristic time 

A k (t) ~ e^-V^o (17) 

where the power of e is due to the locality of the interactions. The behavior ( |TTD will last 
up to the time Tjy at which the perturbation at small scale saturates, i.e. = A sat . Then 
a second regime drived by the variable xjy-i, which has now the fastest exponential growth 
of errors, sets in. When this second regime holds, the incertitudes at scales k < N — 1 
grows as ~ e AiV ~ 1 * until, at time T/v_ 1; the variable xn-i saturates and xn-2 dominates, and 
so on. Such an argument suggests that a perturbation Aj at the largest scale x± follows 
different exponential laws with rates (k = l,---,n) in each different regimes, with a 
global (envelope) evolution which appear very similar to a power law. 

Figure [I] shows the error observed at large scale Aj as a function of time in a numerical 
experiment realized with a chain of N = 10 maps. The coupling function is chosen linear, 
f(x) = x, the coupling parameter is e = 10 -4 , the largest characteristic time is T = 200. 
The initial conditions are random and the initial perturbation on the small scale is 5q = 10~ 8 . 
Figure |I| shows the error observed at large scale Ai as a function of time. Some exponential 
regimes are still recognizable, the straight lines with different slopes in the linear-logarithmic 
plot. 
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The observed behaviour can be well understood by means of a quasi-linear analysis which 
also allows for some analytical estimates. The evolution with (|T5| ) and linear f(x) leads to 
a linear evolution for the perturbation, 



A(t + 1) = AA(t) 



(18) 



where A is the N x N Jacobian matrix 
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h e . 
l 2 e . 



..0 
..0 



, ... l N , 



(19) 



Let us denote with T k the time at which scale k saturates, i.e. A k (T k ) = 1, and with M k the 
period (number of steps) during which scale k dominates the dynamics, i.e. M k = T k — Tk+i- 
If we suppose that the main contribution to the growth of perturbations during period M k 
is given by the faster scale A^, which is correct whenever e is very small, we can write for 
any i < k 



Ai(T k ) = A t (T k+1 ) + (A k ) iyk A k (T k+1 



(20) 



while Aj(Tfc) = A sat for i > k. The matrix A M can be easily computed and at the leading 
order is given by 

3 MA fe _ e MM 



{A M \ k = e 



k-i' 



(21) 



The saturation periods are given by 



Afc A k (T k+ i) 

so that from (gOj) , (gTJ), (|22p and using the fact that e MkXk A k {T k+1 ) ~ A sat we obtain 

Afc e 



(22) 



(23) 



Figure ^ shows the result obtained averaging over 1000 initial perturbations. Scattered 
markers represent the values of Ai(T k ) at the different saturation times, obtained according 
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to the quasi-linear analysis. The agreement with the direct simulation of ([14]) is quite good. 
The global behavior can be fitted by a power law, whose slope can be roughly estimated by 
the following argument. From (Effi), one has that the periods Mk scale, a part logarithmic 
corrections, as Mk ~ A^ 1 ~ a~ k . From ([H]) the growth of Ai during period Mk is estimated 
to be ~ e fe . Then, we expect A 1 (t) ~ t 7 with 7 = — loge/ log a. 

We have also performed a second numerical experiment with a stronger coupling e = 10~ 2 
and quadratic coupling function f(x) = x 2 . The initial error is So = 10~ 3 , while all the other 
parameters in the model are unchanged. We take an average of 1000 runs of 2000 steps in 
order to get a good statistics. The results are reported in figure |3], where one sees that the 
error on large scale Ai(£) has a complex behaviour resulting from the combinations of several 
time-scales, which are now not well separated. However, even in this case, it is possible to 
determine an apparent power-law regime for intermediate times. 

The system of coupled maps with different chaotic characteristic time teaches us that it 
is possible to have non trivial time evolutions of non-infinitesimal perturbations at slowest 
scale, and that these can be fitted by power laws, although they are actually generated by 
saturation processes. In this kind of systems, if we are interested in long time predictions, 
the maximum Lyapunov exponent is of little or no importance and one should consider all 
the time scales present in the system. 



IV. PREDICTABILITY IN A SHELL MODEL OF TURBULENCE 

In this section we discuss the problem of predictability in a model for the energy cascade 
in three dimensional fully developed turbulence fi"5|- |18f| . The model, which is intended to 
mimic the Navier-Stokes equations, is defined in the Fourier space as follows. The Fourier 
space is divided into iV shells labelled by the wavevector modulus k n = k 2 n , where ko is 
an arbitrary constant, each one containing all the wavevector k with modulus ko2 n < \k\ < 
ko2 n+1 . The velocity difference over the length scale ~ k~ x is given by a complex variable 
u n representing the Fourier components of the velocity field. The evolution is given by the 
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set of N ordinary differential equation 
d 

—u n = -vk 2 n u n + ig n + f n (24) 

1 1 

9n — K n U n+l U n+2 — — K n -1 u n-l U n+1 ~ o K n-2 U n _ 2 U n-l 

where f n is an external forcing, acting on large scale, necessary to have a stationary state, 
and v the viscosity. The main differences with the Navier-Stokes equations are: 

• i) the wave-vectors and the velocity fields Un are scalars; 

• ii) there are only nearest and next-nearest neighbour interactions among shells. 

The first point implies that it is not possible to have geometrical structures, all informations 
on phases being lost. The second point does not represent a strong limitation as long 
as the energy cascade in Fourier space is local, with exponentially decreasing interactions 
among shells. This is rather sensible for three-dimensional turbulence, but not for the two- 
dimensional case 

Despite the fact that the time evolution governed by (|2~4])-(|2l)D spends long times around 
an unstable fixed point given by the Kolmogorov law \u n \ oc k^ 1 ^ 3 ||18|| , it exhibits chaotic be- 
haviour on a strange attractor in the 2N dimensional phase space, with maximum Lyapunov 
exponent roughly proportional to v~ l l 2 . The velocity structure functions (\u n \ p ) ~ k~^ p have 
non-linear exponents ( p very similar to those of real turbulence [16| . 



Our study of predictability is based on the comparison of the temporal evolution of pairs 
of different realizations of the velocity field, say u n and v n . Both fields evolve according to 



(p4|)-(|25|) form initial conditions such that: 



i) the energy spectra of u n and v n at the initial time are equal; 

ii) u n and v n at initial time differ only on small scales, corresponding to wavenumber 



kri ^ k n ^ . 



, u n for n < n E 
v n ={ (26) 
e lWn u n for n > ue 
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where w n is a random number uniformly distribited in the range [0, 6\. 

By changing the value of 9 we can modify the correlation between the two field. The 
extreme case 6 = 2it, i.e. completely uncorrelated variables for n > n&, corresponding to 
the Lorenz choice discussed in section ||. 

Previous works |J investigated the growth of different small perturbations in shell mod- 
els. Here we want to study this issue for a finite perturbation at initial time. 

In the numerical simulation we have taken as the reference field u n in (|26D a solution 
of (|24|) -(^5|) obtained from a long simulation starting from a random initial condition. The 
forcing f n was taken constant and equal to 

/n=(l+i)7*n,4 (27) 

with 7 = 0.005. The simulations were done mainly for two different systems: (A) N = 19 
and v = 10~ 6 ; (B) N = 27 and v = 10~ 9 . We used the leap-frog method of integration with 
time step 10~ 3 in the case (A) and 2 x 10~ 5 in the case (B). The average is over 10 3 different 
realizations of the field v n . All the figures will be for the case (A), the case (B) being similar. 

In order to understand whether intermittency has any effect on the qualitative features, 
we also performed a detailed study in terms of closure approximation, where, by construction, 
one considers averaged quantities neglecting intermittency. We thus developed the standard 
eddy damped quasi-normal markovian (EDQNM) approximation |14| for the shell model 
(see appendices [A] and |B|). 

Among the various quantities that can be computed, we focused on the error-spectrum 
at wavevector k n 

^n = \ (K - v n ) « - <)) = (E n - Re W n ) (28) 

where E n = (u n u* n ) = (v n v*) is the energy of the two fields and W n = (u n v*) is the overlap 
energy at scale n. 

In Figs. f| and |5] the error energy spectrum A n is showed at different times for 9 = 2n 
and 10~ 4 , respectively. 
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In figure |] we see that at the beginning we have a very fast growth until A n reaches the 
saturation (i.e. A ra = E n ) at large n. Then one has a sort of inverse cascade on the error. 
This is also the main feature one observes in the case of strong initial perturbation (figure 

D- 

In order to compare the error growth in the shell model with the behaviour observed in 



the toy model of section |T|, we also compute the (global) growth of the moments of the 
difference field, defined as 

W2\ 



(H 9 > 

and the large-scale error 



£2A r 



(29) 



A 6 + A 7 + A 8 



1/2 



(30) 



3 

Figures ^| and [7| show the growth of the moments of (\5u\ q ) 1 ^ q respectively for 6 = 2n 
and 9 = 10 -4 . For the first case we observe at intermediate times a power law behaviour in 
agreement with the Lorenz prediction. The intermittency produces anomalous scaling, i.e. 
(\Su\ q ) ~ t qa{q) where a(q) is a non increasing function of q. Let us note that a(2) = 1/2, 
as in Lorenz prediction where no intermittent effects are taken into account. In the case of 
small initial perturbation (figure one can recognize an initial exponential growth until the 
error saturates at small scales. Then we find a kind of power law behaviour. 

For a qualitative comparison with the results discussed for the coupled maps, we show 
in figure [8] the growth of 5Uq. We can recognize an initial exponential growth driven by the 
smallest scale; in later times all the time scales interplay leading to a non simple function 
which can be misinterpreted as a power law. 

The closure approximation for the shell model leads to a qualitaive similar behaviour. 
The basic idea of closure is quite straightforward: one writes a Reynolds hierarchy for the 
moments of the shell variables, and truncate the chain of equations at the lowest sensible 
order. The full derivation is given below in appendices [A] and |B|. The closure equation thus 
gives information of the behaviour of E n and W n as functions of time. The results are similar 
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to those obtained by the direct integration (and ensamble average) of the shell model, as 
is evident by the direct comparison of figure |] with figure ^ and figure [10] with |5]. This is 
a clear evidence that the relevant mechanism, at least at a qualitative level, is due to the 
existence of many characteristic times and not to intermittency effects. 

V. DISCUSSION 

The predictability problem can be formulated in terms of the Lyapunov exponent, or 
the effective Lyapunov exponent, only for infinitesimal perturbations. On the other hand 
in systems with many different characteristic times and finite perturbations the non-linear 
effects are very relevant. Direct consequence is that the growth of the perturbation is, a 
part at very short times, not exponential but rather similar to a power law. Actually it is 
a sequence of exponentials with varying rates. This mechanism has an important practical 
relevance: the predictability time for finite perturbations can be much larger than the inverse 
of the Lyapunov exponent ([[]). 

For example, in three dimensional turbulence for a "physical" perturbation localized 
at the initial time at the Kolmogorov lenght i], Sv = 0(v(rj)), the predictability time is 
T ~ L/V, independent of the Reynolds number. On the contrary, the predictability time 
for an infinitesimal perturbation is proportional to the inverse of the maximum Lyapunov 
exponent A, which grows as a power of the Reynolds number. Hence the predictability time 
decreases for increasing Re. 

These results give a strong indication that despite the presence of strong chaos, realistic 
situations have relatively long predictability times. 

In conclusion we have showed that the concepts of Lyapunov exponents, effective Lya- 
punov exponent, Kolmogorov entropy give insufficient informations on the chaotic behaviour 
of extended systems, and it is urgent the identification of new indicators. 
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APPENDIX A: 



Here we derive the equations for the energy of the field in the eddy damped quasi-normal 



markovian approximation (EDQNM) for the shell model. For more details see e.g. [14,19|. 

We are interested into the energy of the shell n given by E n = (u n w* ) . Differentiation 
with respect to time of E n , and use of ( p^) - (|25|) , leads to the evolution equation for E n for 
the shell model: 

^E n = -2 v k 2 n E n + i(g n u* n ) - i(u n g* n ) + 2 e 5 nA . (Al) 
In order to have a constant energy input e we assumed a forcing 



fn=j\e5 nA (A2) 



in (|24l). This equation is not closed since the third and fourth terms on r.h.s involve third 
order correlation functions. For example 

(g n u* n ) = M<<+i<+2> - 2 ^-i«-i<<+i) -^i-2«-2<-i<) (A3) 

and similarly the other. The third order cumulants obey a differential equation involving 
fourth order correlation functions, e.g. 

4 (u n u n+1 u n+2 ) = - is(k 2 n + k 2 n+1 + h 2 n +2 ) (u n u n+1 u n+2 ) 

+i{g n u n+ iu n+2 } +i(u n g n+ iu n+2 ) + i(u n u n+1 g n+2 ), 
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and so on. Thus in order to have a closed problem we have to truncate the hierarchy. This 
is done assuming a quasi-normal probability distribution for the field, and factorizing the 
fourth order correlation functions as 



(xiX 2 X 3 X4) = (X1X2) (X3X4) + (X1X3) (X2X4) + (xiXi) (x 3 x 2 ). 



(AS) 



By means of (|A5|) , and the assumption of isotropic turbulence 



( u n <) = S nm E n , (u n u m ) = 



(A6) 



it is easy to evaluate the terms involving the fourth order correlation functions in the equa- 
tions for the third order ones, and close the hierarchy. For example for Eq. (|A4j) we have, 



Inserting 



\9n u n+l u n+2 ) — k n E n+ i E n+2 

(u n 9n+l Mn+2) — ~\k n E n E n+2 

(u n u n+ i g n +2) = —\k n E n E n+ i. 

into ( |A4| ) we obtain 



(A7) 



37 + ^(^n + kn+1 + kn+2) 



(u n u n+1 u n+2 ) = iS(n,t) 



(A8) 



with 



1 1 

S{n, t) = k n E n+ i E n+ 2 — — k n E n E n+2 ~ — k n E n E ri+ \. 



(A9) 



Similar equations hold for the third order cumulants entering into Eq. (|A3|). 

The quasi-normal approximation (|A8|) it is known to give unphysical results [ 19[| . In 
particular it does not produce positive definite energy spectra. The problem has been cured 



by introducing the so-called "eddy damped" approximation. One replaces (|A8|) with 

+ U (kn + ^n+l + kn+2) 

+ /in + /i„+i + /i„+2 (u n u n+1 u n+2 } =iS(n,t) 

where 



(A10) 



= KK, E n) =ak„ El /2 . 



(All) 
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Dimensionally k n E^ 2 is an inverse time, the turn-over time at shell n. We have one free 
parameter, the dimensionless constant a. It should be adjusted such that the spectrum is 
as similar as possible to the spectrum obtained in simulations of the full equation. We will 
return to this point later. 

Equation ( |A10|) can be easily integrate, and the solution reads 



(U n U n + 1 U n+2 ) = I f dt' e -K^+^ + 1 +^ +2 )+Mn+Mn + l+^ +2 ] (t-f) gfa (A12) 

Jo 



It has been proved that if we assume that S(n, t) does not vary significantly in the range 
where the exponential in ( |A12| ) is substantially different from zero, then this is a suffi- 



cient condition for the positiveness of the energy spectra ]nj. This assumption is called 
Markovianization. Therefore the third order correlation function at time t in the EDQNM 
approximation reads 



(u n u n+ i u n+2 ) = iO(n, t) S (n, t) (A13) 



where 



Jo 

I _ e -H k l + k l + l+ k l + 2 )+t 1 ™+Hn + l+Hn + 2}t 



(A14) 



v{k 2 n + k 2 n+l + k 2 n+2 ) + /i n + fi n+1 + fi n+2 ' 
Similar equations hold for the other third order correlation functions. 

We can now go back to ( |A1| ) and write down the equation for the energy E n in the 
EDQNM approximation 



(i + 2uk 2 n ) E n = 2 [k 2 n e(n,t)S(n,t) 



kl^ein-l^Sin-^t) (A15) 



.1 U2 

2 



■±k*_ 2 9(n-2,t)S(n-2,t) +2e6 nA . 



The quasi-normal ansatz implies the absence of intermittency corrections. This is an 
essential limitation of all closure theories. The energy spectrum of the shell model in the 
EDQNM approximation must therefore obey E n ~ C(a)e 2 ^ 3 k~ 2 ^ 3 in the inertial range. The 
undetermined function C(a) is the Kolmogorov constant. 



On the other hand it has become clear in several independent investigations that in- 
termittency corrections exist in shell models. The energy spectrum is therefore in reality 
more closely described by E n ~ F(e)k~^ 2 , where the exponent (2 has been estimated to be 
0.70 ||16| . The function F that gives the prefactor to the power law in the inertial range 
should not depend on viscosity, but depends on the forcing through e, the mean dissipation 
of energy per unit time, or, equivalently, the mean energy input into the system from the 
force. In a really large inertial range the two power-laws are not good approximations to 
one other. The best that can be done is to demand that the spectra agree as closely as 
possible at the upper end of the inertial range. The disagreement at the lower end will then 
be approximately (&l//k)) < ' 2 ~ 2 ^ 3 - This is not a very large discrepancy. Assuming an inertial 
range of twenty shells, that is a scale separation of 10 6 , and (2 — 0.70, the mismatch is only 
a factor 1.5. In practice a number of a's have been tried and the two spectra compared until 
a reasonable agreement is achieved. For a = ??? we obtain C(a) = 1.5 which is the value 
observed both in simulations of the shell model and in experiments |TL9 . 



APPENDIX B: 

Here we derive the equations for the energy of the field difference in the eddy damped 
quasi-normal markovian approximation (EDQNM) for the shell model. The procedure is 
similar to that described in App. [A] for the energy, so we only report the main equations. 

We consider to independent realizations of the field, u n and v n , with the same energy 
spectrum E n , both evolving according to (^) and ( flop . For simplicity of notation, the 
equation of motion of the field v n is rewritten as 

A Vn = - v k 2 n + ih n + /„ 

h n — k n V * +1 V^ + 2 ~ 2 ^n-l v n-l V n+1 ~~ 2 V n-1 (^1) 

Jn — 1 T2" e °n,4 

\ V n\ 

We are interested into the energy of the field difference at the shell n 
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A n = i ((«„ - v n ) « - <)) = (E n - Re W n ) (B2) 

where W n = (u n u*). The evolution equation of W n is easily derived by differentiation with 
respect to time and reads, 

^-W n = -2uk 2 n W n + i{g n v*) - i{u n h* n ) + 2e 5 nA (B3) 
at 

As the case of E n this equation is not closed since it involves third order cumulants of u n 
and v n . For example 

(9nV*) = k n (v* n U* n+l ul +2 ) - 2*Vi-l«-l<«Jl+l) - ^n-2«_ 2 <-lO ( B4 ) 

and similarly the other. As done for the energy E n , the hierarchy is closed the fourth order 
correlation functions. With a calculation similar to that of App. |A], one obtains 
d 



ft + v{K + k 2 n+1 + k 2 n+2 ) 



( V n U n+1 U n+2) ~ 

(B5) 

-ik n W* +1 W* +2 + i\ k n W n E n+2 + %\ k n W n E n+1 



and similar equations for the others. 

We then perform the eddy damped and markovian approximation, i.e. add a damping 
term (/i n + fi n +i + ^n+2)W n to the l.h.s, where \i n is given by (|A11|) , and we integrate the 
resulting equation in the markovian approximation - see App. [A] for more details. From 
(|B5D, it follows 



« <+i <+ 2 ) = -W(n, t) k n [k n W: +1 W* +2 

(B6) 

-\ k n W n E n+2 -\k n W n E n+ l] 



with 9(n,t) given by ( |A14j ). The other cumulants leads to similar equations. Collecting all 



the terms one finally obtains the equation for W n in the EDQNM approximation 

(i + 2vk 2 n ) W n = 2 [k 2 n 9(n, t) (W: +1 W* n+2 - \ W n E n+2 - \ W n E n+1 ) 

-\ kl_ x 6(n - 1, t) (W n E n+l - \ W*_ x W r * l+1 - \ E n -, W n ) 

(B7j 

_ 1 k 2 n _ 2 6(71 - 2, t) (E n -! W n - \ E n . 2 W n - \ W* n _ 2 W*_ x ) 
+2eS n A. 
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FIGURES 

FIG. 1. Error growth Ai as a function of t for e = 10 -4 . The straight segments correspond to 
exponential growth driven by different dominant scales. 

FIG. 2. (Ai) as a function of t for e = 10~ 4 . The average is taken over 1000 initial configura- 
tions. The symbols (dyamonds) are the estimated growth of the largest scale Ai during different 
periods according to the quasi-linear approximation. Dashed line: power law expected by the 
simple argument described in the text. 

FIG. 3. Error growth Ai for a non-linearly coupled system with e = 10~ 2 . 

FIG. 4. Error spectrum A n as a function of n at different time for case (A) and 9 = 2ir. The 
lines are taken at times:0, 1, 2, 3, 4, 5, 10, 15, 20, 25 and 30 seconds after the perturbation has 
been done. The dashed line is the spectrum E n . 

FIG. 5. Error spectrum A n as a function of n at different time for case (A) and 9 = 10~ 4 . The 
lines are taken at times:0, 5, 10, 15, 20, 25, 30, 35, 40, 45 and 50 seconds after the perturbation 
has been done. The dashed line is the spectrum E n . 

FIG. 6. Moments of the difference field (\8u\*) 1 l« (q = 1,2,3,4) as a function of t for case (A) 
and 9 = 2tt. 

FIG. 7. Moments of the difference field (\5u\ q ) 1 / q (q = 1, 2, 3,4) as a function of t for case (A) 
and 9 = 10~ 4 . 

FIG. 8. SUq as a function of t for the case (A). 

FIG. 9. Error spectrum A n as a function of n at different time for the EDQNM approximation 
and 9 = 10~ 4 . The lines are taken at time interval 0.02 seconds from the the perturbation, first 
line. The dashed line is the spectrum E n . 
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FIG. 10. (|<5u| 2 ) as a function of n at different time for the EDQNM approximation and 
6 = 10" 4 . 
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